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Abstract 

<N ■ 

J> ■ We generalize previous work on early photoionization to CDM 

^] . models with Q < 1. Such models have received recent interest because 

the excess power in the large-scale galaxy distribution is phenomeno- 
' logically fit if the "shape parameter" T = hQ® ~ 0.25. It has been 

argued that such models may require early reionization to suppress 
\ small-scale anisotropics in order to be consistent with experimental 

| data. We find that if the cosmological constant A = 0, the extent of 

O this suppression is quite insensitive to 0,q. Given a erg-normalization 

q . today, the loss of small-scale power associated with a lower Qq is par- 

S-) | tially canceled by higher optical depth from longer lookback times and 

by structures forming at higher redshifts before the universe becomes 
curvature-dominated. The maximum angular scale on which fluctu- 
. ^ | ations are suppressed decreases when is lowered, but this effect is 

^ ' also rather weak and unlikely to be measurable in the near future. For 

flat models, on the other hand, where Ao = 1 — ^o> the negative effects 
of lowering dominate, and early reionization is not likely to play a 
significant role if <C 1. The same applies to CDM models where the 
effective T is lowered by increasing the number of relativistic particle 
species. 
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1 Introduction 



Early reionization of the universe would affect fluctuations in the cosmic mi- 
crowave background radiation (CBR) on degree-scales and sub-degree scales. 
Considerable effort is now being devoted to measuring such fluctuations, so 
a detailed knowledge of the ionization history is desirable in order to give the 
experimental results their proper cosmological interpretation. 

It is generally believed that the intergalactic medium (IGM) was reion- 
ized at some time in the past. The main reason for this is the absence of 
a Gunn-Peterson trough in the spectra of high redshift quasars (Gunn & 
Peterson 1965; Steidel & Sargent 1987; Webb et al. 1992), indicating that 
reionization occurred at least as early as z = 4. For CBR applications, the 
relevant question is not whether reionization occurred, but when it occurred. 
For instance, to provide an optical depth to scattering exceeding 20%, which 
suffices to reconcile CDM with all observational CBR limits (Sugiyama, Silk 
& Vittorio 1993), reionization must have occurred by redshift z = 30, adopt- 
ing the nucleosynthesis value of the baryon density. 

In Tegmark, Silk & Blanchard 1994 ("Paper 1"), early photoionization 
of the intergalactic medium was discussed in a fairly model-independent 
way, in order to investigate whether early structures corresponding to rare 
Gaussian peaks in a CDM model could photoionize the intergalactic medium 
sufficiently early to appreciably smooth out the microwave background fluc- 
tuations. In this paper, the results of Paper 1 will be generalized to Q < 1 
models with non-zero cosmological constant A. Essentially all the notation 
used in this paper was defined in Paper 1 and, in the interest of brevity, some 
of the definitions will not be repeated here. 

Just as in Paper 1, our basic picture is the following: An ever larger 
fraction f s of the baryons in the universe falls into nonlinear structures and 
forms galaxies. A certain fraction of these baryons form stars or quasars 
which emit ultraviolet radiation. Some of this radiation escapes into the 
ambient intergalactic medium (IGM), which is consequently photoionized 
and heated. Due to cooling losses and recombinations, the net number of 
ionizations per UV photon, f ion , is generally less than unity. 

The results that we present here generalize the previous work to the case 
of an open universe. Lowering fl has four distinct effects: 

1. Density fluctuations gradually stop growing once z ^ f2g 1 . Thus given 
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the observed power spectrum today, a lower Q implies that the first 
structures formed earlier. 

2. Matter-radiation equality occurs later, which shifts the turning-point 
of the CDM power spectrum toward larger scales. This means less 
power on very small scales (such as ~ 10 6 M Q ) relative to the scales at 
which we normalize the power spectrum (namely galaxy cluster scales, 
~ 8/i _1 Mpc or even the much larger COBE scale). One consequence is 
that the first structures form later. 

3. The lookback time to a given ionization redshift becomes larger, result- 
ing in a higher optical depth. 

4. The horizon at a given ionization redshift subtends a smaller angle on 
the sky, thus lowering the angular scale below which CBR fluctuations 
are suppressed. 

Thus in terms of the virialization redshift z V i r defined in Paper 1, the redshift 
at which typical structures to go nonlinear, effect 1 increases z V i r whereas ef- 
fect 2 decreases z V i r . These two effects influence f s , the fraction of baryons in 
nonlinear structures, in opposite directions. As to effect 2, it should be noted 
that this applies not only to CDM, but to any spectrum that "turns over" 
somewhere between the very smallest nonlinear scales (~ 1O 6 M ) and the 
very largest (~ 1O 21 M ) scales at which we COBE-normalize. Yet another 
effect of lowering Q is that the ionization efficiency /j ora drops slightly, at 
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most by a factor f2 . This is completely negligible compared to the above- 
mentioned effects, as z ion depends only logarithmically on the efficiency. 

In the following sections, we will discuss each of these four effects in 
greater detail, and then compute their combined modification of the ioniza- 
tion history in a few scenarios. 



2 The Boost Factor 

When curvature and vacuum density are negligible, sub-horizon-sized density 
fluctuations simply grow as the scale factor a oc (1+z)^ 1 . Thus at early times 
z ^> Qq 1 , we can write 

o = — — o 

1 + z 
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for some function B independent of z that we will refer to as the boost factor. 
Clearly 5(1,0) = 1. Thus if certain structures are assumed to form when 8 
equals some fixed value, then given the observed power spectrum today, the 
boost factor tells us how much earlier these structures would form than they 
would in a standard flat universe. The boost factor is simply the inverse of 
the so called growth factor, and can be computed analytically for a number 
of special cases (see e.g. Peebles 1980). For the most general case, the fit 



^ /7 -Ao+(l + ^ 




is accurate to within a few percent for all parameter values of cosmological 
interest (Carroll et al. 1992). The exact results are plotted in Figure [I] for 
the case Ao = and the flat case Ao = 1 — fio- Since we will limit ourselves 
to these two cases, the simple power-law fits 
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r B(n ,o) «a 

is(fio,l-fio) ~^c7°- 21 , 

which are accurate to within 1% for 0.2 < Qq < 1, will suffice for our pur- 
poses. Note that the standard rule of thumb that perturbations stop growing 
at 1 + z ~ , indicating B cx VLq , is not particularly accurate in this con- 
text. 



3 The Power Spectrum Shift 

As mentioned above, lowering M7 causes the first structures go nonlinear at 
a later redshift. This is quantified in the present section. 

The standard CDM model with power-law initial fluctuations propor- 
tional to k n predicts a power spectrum that is well fitted by (Bond & Efs- 
tathiou 1984; Efstathiou et al. 1992) 



(l + fag+^^ + ^g) 2 ] 1 " 13 ) 271 ' 13 ' 

where a = 6.4, b = 3.0, c = 1.7, q = (l/i -1 Mpc)A;/r and the "shape pa- 
rameter" T will be discussed below. Although this fit breaks down for scales 
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comparable to the curvature scale r curv = -f^o" 1 ]! — ^ | _1 , it is quite accu- 
rate for the much smaller scales that will be considered in the present paper. 
Rather, its main limitation is that it breaks down if Q is so low that the 
baryon density becomes comparable to the density of cold dark matter. Thus 
for flj, pa 0.05, the results cannot be taken too seriously for Q < 0.2. We will 
limit ourselves to the standard n — 1 model here, as the tilted (n < 1) case 
was treated in Paper 1 and was seen to be be essentially unable to reionize 
the universe early enough to be relevant to CBR anisotropics. The same 
applies to models with mixed hot and cold dark matter. 
Let us define the amplitude ratio 

fl( r, ri ,r 2 )^, 

a{r 2 ) 

where cr(ri) and a(r2) are the r.m.s. mass fluctuation amplitudes in spheres 
of radii r\ and r 2 , i.e. 



POO 

a(r) 2 oc / P(k) 
Jo 



sin kr cos kr 1 " 



(kr) 3 (kr)' 



dk. 



As in Paper 1, we normalize the power spectrum so that a(8h~ l Mpc) equals 
some constant denoted by as, and M c will denote the characteristic mass 
of the first galaxies to form. The corresponding comoving length scale r c is 
given by M c = ^nr^p, where p is the mean density of the universe. Thus 
given as, what is relevant for determining when the first galaxies form is the 
amplitude ratio 

R(T,r c ,8h- l Mpc). 

This ratio is computed numerically, and the results are plotted as a function 
of r in Figure H for a few different values of the cutoff mass M c . It is easy to 
see why the amplitude ratio increases with T, since on a logarithmic scale, a 
decrease in T simply shifts the entire power spectrum towards lower k, thus 
decreasing the amount of power on very small scales relative to that on large 
scales. The fit 

R{F, r c , 8/i _1 Mpc) ^3 + 7.1 ln(l/i -1 Mpc/r c )r 
is accurate to within 10% for 0.05 < T < 2 and lOOpc < r c < lOOkpc. 
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4 The Optical Depth 



Since a lower Q implies a larger \dt/dz\ and an older universe, the optical 
depth out to a given ionization redshift Zi on is greater for small Qq. For a 
given ionization history x{ z )i the optical depth for Thomson scattering is 
given by 

t(z) =t*L z , xQW, 



3^6 
8vr 



1 - (l - *g) he] - 0.057Ml b , 



where we have taken the mass fraction of helium to be fn e ~ 24% and 
assumed \He ~ X-, i- e - th & t helium never becomes doubly ionized and that 
the fraction that is singly ionized equals the fraction of hydrogen that is 
ionized. The latter is a very crude approximation, but has the advantage 
that the error can never exceed 6%. If the universe is fully ionized for all 
redshifts below z, the integral can be done analytically for A = 0: 



2t* 



2 - 3fi + (Poz + 3fi - 2)yJl + Q z 



,MV 3 / 2 



0.038^^- (1) 



for z ^> Qq 1 . As is evident from the asymptotic behavior of the integrand, r 
is independent of Ao in the high redshift limit. Thus optical depth of unity 
is attained if reionization occurs at 



-2/3 



1/3 



Z « 92 fin 

V0.03; 

Besides fluctuation suppression, polarization of the CBR is another in- 
teresting probe of reionization (Crittenden, Davis & Steinhardt 1994). Also 
for this application, the optical depth is a key parameter. 



5 The Angular Scale 

It is well known that reionization suppresses CBR fluctuations only on angu- 
lar scales below the horizon scale at last scattering. Combining the standard 
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expressions for horizon radius (e.g. Kolb & Turner 1990) and angular size, 
this angle is given by 



9 = 2 tan 



-i 



sqrtl + zn 3 J 2 /2 



n oZ - (2 - n )(y/i + n z - i) 



(2) 



For z ^> Q , this reduces to 




(3) 



but as is evident from Figure |3], this is quite a bad approximation except for 
z ^> 100. If we substitute it into equation ([]]) nonetheless, to get a rough 
estimate, we conclude that optical depth unity is obtained at an epoch whose 
horizon scale subtends the angle 



9 « 12° 



6"0 



1/3 



0.03 



(4) 



i.e., the dependence on all three of these cosmological parameters is relatively 
weak. As discussed in Paper 1, fluctuations on angular scales much smaller 
than this are suppressed by a factor 



P(z) 



the opacity, which is the probability that a photon was Thompson scattered 
after redshift z. Its derivative, the visibility function f z = dP/dz, is the 
probability distribution for the redshift at which last scattering occurred, 
the profile of the last scattering surface. The angular visibility function 



fe(0) 



dP« 



de 



de 



dz 



dPs 
dz 



is plotted in Figure |] for the case where the universe never recombines (the 
curves for the more general case with reionization at some redshift Zi on can be 
read off from Figure |] as described in Paper 1). These functions give a good 
idea of the range of angular scales on which suppression starts to become 
important. In plotting these curves, the exact expression (||) has been used, 
rather than the approximation @. It is seen that the qualitative behavior 
indicated by equation (|j) is correct: as Qq is lowered, the peak shifts down 
toward smaller angular scales, but the f2 -dependence is quite weak. 
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6 Cosmological Consequences 



We will now compare the effect of lowering T in three cosmological models. 
The first model, which will be referred to as "open CDM" for short, has 
Ao = 0. The second model, referred to as "ACDM", has Ao = 1 — £V The 
shape parameter essentially tells us how early the epoch of matter-radiation 
equality occurred, and is given by 



hQr 



9* ^ 1/2 



3.36 



where g* = 3.36 corresponds to the standard model with no other relativistic 
degrees of freedom than photons and three massless neutrino species. In 
open CDM and ACDM, we have g m = 3.36, so that Qq = T/h. The third 
model, referred to as rCDM (Dodelson, Gyuk & Turner 1994), has A = 
and Q = 1, and achieves a lower value of T by increasing instead. 

Including the effect of the boost factor, equation (9) in Paper 1 becomes 

1 + z lon = R(m , r c , S/i^Mpc) B(Q , A ) erfc" 1 

Or 



fuvppfion _ 



The ionization redshift Zj on is plotted as a function of the shape parameter 
in Figure |5| for the various scenarios specified in Table |l| It is seen that for 
the open model, the dependence on Qq = T/h is typically much weaker than 
the dependence on other parameters. One reason for this is that changes 
in the boost factor and the amplitude ratio partially cancel each other. For 
ACDM, the ^-dependence is stronger, since the boost factor is weaker. In 
the rCDM model, the dependence on T/h is even stronger, as there is no 
boost factor whatsoever to offset the change in the amplitude ratio. 

The scenarios in Table [I] are similar to those in Paper 1. In the one labeled 
"very optimistic" , the high value for f uvpp , the net number of produced UV 
photons per proton, is obtained by assuming that the main source of ionizing 
radiation is black hole accretion rather than conventional stars. Note that 
this speculative assumption still only increases z ion by 3.55/3.00 — 1 ~ 18%, 
the efficiency dependence being merely logarithmic. 

Figure |5], in a sense the most important plot in this paper, shows the 
opacity as a function of T/h for the various scenarios. Because of the increase 
in optical depth due to larger lookback times, the open model now gives 
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Pess. 


Mid. 


Opt. 


Very opt. 


^8 


0.5 


1 


1.1 


1.2 




2.00 


1.69 


1.44 


1.33 


h 


0.5 


0.5 


0.8 


0.8 


M C [M & ] 


10 s 


10 6 


10 5 


10 5 


f f ion f uvpp 


1 


120 


23,000 


10 6 


erfc- 1 [l/2/] 


0.48 


2.03 


3.00 


3.55 


h 2 n b 


0.010 


0.013 


0.015 


0.020 



Table 1: Parameters used 



slightly larger opacities for lower Qq = T/h. However, this dependence is 
seen to be quite week. For ACDM, where the boost factor contributes less, 
the net result is seen to be the opposite; a slight decrease in the opacity 
for lower Qq — T/h. For the rCDM model, where there is neither a boost 
factor nor an increase in the lookback time, this drop in opacity is seen to 
be much sharper. Note that the dependence on other uncertain parameters, 
summarized by the four scenarios in Table [1[ is quite strong. Indeed, this 
dependence is stronger than the effect of moderate changes in T/h, so in the 
near future, it appears unlikely that opacity limits will be able to constrain 
the shape parameter except perhaps in the rCDM model. 

The rCDM situation is summarized in Figure |7]. To attain at least 50% 
opacity, hQb must lie above the heavy curve corresponding to the scenario 
in question. On the other hand, nucleosynthesis (Smith et al. 1993; Walker 
et al. 1991) places a strict upper bound on this quantity if we assume that 
h > 0.50 It is seen that a shape parameter as low as T ~ 0.25, which would 
match large-scale structure observations (Peacock & Dodds 1994), is quite 
difficult to reconcile with these two constraints. 

7 Discussion 

Lowering T is an attractive resolution of the problem that arises in reconcil- 
ing the observed structure in the universe on large scales with observations 

2 The rCDM model also alters the nucleosynthesis process (Dodelson, Gyuk & Turner 
1993), but this can only marginally relax the bounds unless the r neutrino is in a mass 
range incompatible with rCDM (Gyuk & Turner 1994). 
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on megaparsec scales. The empirical power spectrum is well fit by T 0.25 
(Peacock & Dodds 1994). Kamionkowski & Spergel (1993) have found that 
primordial adiabatic fluctuations in an open universe with Q 0.3 are rec- 
oncilable with large-scale CMB anisotropy. On degree scales Kamionkowski, 
Spergel & Sugiyama (1994) require reionization with optical depth r ~ 1 in 
order to reconcile the low density open model with recent experimental lim- 
its, if the lowest limits are adopted. In a low Q ACDM model, the situation 
is not so critical, but reionization is required if the lowest limits (SP91) are 
adopted on degree scales (Gaier et al. 1992); r ~ 0.5 suffices however. A 
similar but slightly more favorable situation occurs in a rCDM model, where 
T = M)o(<7*/3.36) -1 / 2 is reduced by increasing g* by a factor of ~ 4, but 
some reionization is still required to match SP91. 

We have found that reionization giving r in the range 0.5 to 1 is readily 
produced and even natural in open models. This is because of the early 
formation of structure in combination with the increased age of the universe, 
effects which compensate for the flattening of the power spectrum due to 
the delay in matter domination. However, the ACDM and rCDM models 
with low Q fare less well in this regard, since the loss of small-scale power 
is not balanced by significantly earlier structure formation. With in the 
range given by standard nucleosynthesis, a significant optical depth r > 0.5 
is difficult to attain in either the ACDM or rCDM scenarios. 

The authors would like to thank Wayne Hu, Lloyd Knox, Bernard Sadoulet 
and Douglas Scott for many useful comments. This research has been sup- 
ported in part by a grant from the NSF. 
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0.2 0.4 0.6 0.8 1 



Figure 1: The boost factor. 
The boost factor B(£l , A ) is plotted as a function of f2 for two classes 
of cosmologies. The upper curve corresponds to a standard open universe, 
i.e. Ao = 0, whereas the lower curve corresponds to flat universes with 
Aq = 1 — ^o- 



11 









0.5 1 1.5 

r 



2 



Figure 2: The amplitude ratio. 
The ratio of the fluctuation amplitude on the small scale r c to that at 
8/i _1 Mpc is plotted as a function of the shape parameter T. From top to bot- 
tom, the four curves correspond to scales of 3.5/i _1 kpc, 7.5/i _1 kpc, 16/i _1 kpc 
and 35/i _1 kpc, respectively. For h = 0.5 and Q = 1, these four length scales 
correspond to the masses 1O 5 M , 1O 6 M , 1O 7 M and 10 8 M Q . The weak ad- 
ditional dependence on Q /h that would result from holding M c rather than 
r c fixed is clearly negligible, since as can be seen, M c must vary by an entire 
order of magnitude to offset a mere 20% change in T. 
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Figure 3: The horizon angle. 
The angle in the sky subtended by a horizon volume at redshift z is plotted 
as a function of Qq for the case with no cosmological constant. The solid lines 
are the exact results for the four redshifts indicated, and the dashed lines are 
the corresponding fits using the simplistic approximation 9 « (Hq/z) 1 ^ 2 . 



13 




Degrees subtended by horizon at last scattering 



Figure 4: Visibility functions. 
The angular visibility function for a fully ionized universe is plotted for dif- 
ferent values of ffo and diffuse baryon content hVL b - The left group of curves 
corresponds to hQ b = 0.03 and the right to hQ b = 0.1. Within each group, 
from left to right starting at the lowest peak, f2 =0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 
0.7, 0.8, 0.9 and 1.0. 
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0.2 0.4 0.6 0.8 1 

r/h (=n ) 

Figure 5: The ionization redshift. 
The ionization redshift is plotted for the four scenarios described in Table p], 
the heavy lines corresponding to the optimistic scenario. The solid lines are 
for the open case where A = 0. The dashed lines correspond to the flat case 
where Ao = 1 — Qq. The dotted lines refer to the rCDM model, where Ao = 
and flo = 1. Note that the combination T/h is really an /i-independent 
quantity: for the open and flat cases, it is simply equal to Q Q , and for the 
rCDM case it depends only on the number of relativistic particle species. 
The vertical shaded region corresponds to values of T preferred by power 
spectrum measurements, 0.2 < T < 0.3, when h = 0.5. 
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0.2 0.4 0.6 0.8 1 



r/h (=n ) 

Figure 6: The opacity. 
The opacity, the probability that a CBR photon is Thomson scattered at 
least once since the standard recombination epoch, is plotted for the four 
scenarios described in Table 1. The solid lines correspond to the open case 
where A = 0. The dashed lines are for the flat case where A = 1 — fi . The 
dotted lines correspond to the rCDM model, where Ao = and Q = 1. Just 
as in Figure 5, note that the combination T/h is really an /^-independent 
quantity: for the open and flat cases, it is simply equal to f2 , and for the 
rCDM case it depends only on 7, the number of relativistic particle species. 
The vertical shaded region corresponds to values of T preferred by power 
spectrum measurements, 0.2 < T < 0.3, when h = 0.5. 
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0.2 0.4 0.6 0.8 1 

r 

Figure 7: Reionization in rCDM. 
The two curves show the baryon density required for 50% opacity in rCDM, 
i.e. for reionization to rescatter 50% of the CBR photons. The upper and 
lower heavy curves correspond to the middle-of-the road and optimistic sce- 
narios, respectively. Thus even in the optimistic scenario, 50% opacity can- 
not be obtained outside of the fine-hatched region. The horizontal shaded 
region corresponds to the values of hflf, allowed by standard nucleosynthesis 
(0.01 < h 2 Qt, < 0.015) in conjunction with the constraint 0.5 < h < 0.8. 
The vertical shaded region corresponds to values of the "shape parameter" 
T preferred by power spectrum measurements. Note that these three regions 
do not intersect. 
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